Lorentz Gas at Positive Temperature 
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We investigate the evolution of a particle in a Lorentz gas where the background scatters move 
and collide with each other. As in the standard Lorentz gas, we assume that the particle is negligibly 
light in comparison with scatters. We show that the average particle speed grows in time as f-^/(*+-'') 
in three dimensions if the particle-scatter potential diverges as r~^ in the small separation limit. 
The particle displacement exhibits a universal growth, linear in time and the average speed of the 
atoms. Surprisingly, the asymptotic growth is independent on the gas density and the particle- 
atom interaction. The velocity and position distributions approach universal scaling forms which 
are non-Gaussian. We determine the velocity distribution in arbitrary dimension and for arbitrary 
interaction exponent A. For the hard-sphere particle-atom interaction, we compute the position 
distribution and the joint velocity-position distribution. 

PACS numbers: 05.20.Dd: Kinetic theory, 45.50.Tn: Collisions, 05.60.-k: Transport processes 



The Boltzmann equation [T] is the basic tool in elu- 
cidating the properties of transport phenomena. The 
non-linear integro-difFerential Boltzmann equation is so 
formidable, however, that apart from the equilibrium 
Maxwell-Boltzmann distribution [2] there are essentially 
no solutions to the Boltzmann equation [3]. The stan- 
dard Lorentz gas model where a point particle is elas- 
tically scattered by immobile hard spheres is described 
by the Lorentz-Boltzmann equation U which is linear 
and, not surprisingly, amenable to analytical treatments. 
The Lorentz gas has played an outstanding role in con- 
crete calculations (e.g. of the diffusion coefficient) and 
in the conceptual development of kinetic theory [51 [5]. 
Yet the very applicability of the Boltzmann framework 
to the Lorentz gas is questionable — when the scatters 
are fixed, the molecular chaos assumption underlying the 
Boltzmann equation is hard to justify [SH^. 

If, however, the background particles (atoms for short) 
move and collide with each other, the molecular chaos as- 
sumption holds in the dilute limit and the (properly gen- 
eralized) Lorentz-Boltzmann equation must be applicable 
as long as the mass of the point particle is infinitesimally 
small so that it does not affect the motion of atoms. Sur- 
prisingly this model has not been studied in the context 
of kinetic theory until very recently [9 and many aspect 
of it are still unclear (e.g. density profile). This problem 
is also reminiscent of the model proposed by Fermi [10 to 
explain the acceleration of interstellar particles which has 
been mostly studied using methods of dynamical systems 
(see e.g. [TT] and references therein). 

The emerging behavior of the particle in our model is 
drastically different from the case of the standard Lorentz 
gas with immobile atoms. The average particle velocity 
and and displacement grow linearly in time. Further, af- 
ter a re-scaling with respect to the average quantities, the 
velocity and displacement distributions approach univer- 
sal forms which are not Gaussian. 

The ratio of masses determines the particle equilibrium 
velocity. When this ratio is zero the equilibrium parti- 
cle's velocity is infinite and on the quest to equilibration 
the particle speed increase indefinitely in time. Note that 



in this case the particle carries no kinetic energy and no 
momentum so the conservation of energy and momentum 
do not constrain the particle's velocity. A light particle 
will eventually thermalize with the background atoms at 
some finite velocity and the velocity and displacement 
distribution will become Gaussian. Taking the limit of a 
negligible light particle allow us to push the equilibration 
time to infinity and to observe interesting not equilibrium 
behaviors. Our calculation correctly reproduces the be- 
havior of systems with a mass ratio not strictly zero up to 
some characteristic time where the onset of equilibration 
appears [12]. 

Let us first analyze the velocity distribution. Sup- 
pose the atoms are hard spheres of radius a. The par- 
ticle velocity distribution /(v, t) satisfies the Lorentz- 
Boltzmann equation 



dnP{n)ga''-' / De[/(v',t) - /(v,t)] . (1) 
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dt 



Here P(u) = p(27rr)-'^/2g-«V2T jg ^-^^ Maxwell- 
Boltzmann velocity distribution of the background gas 
(d is the spatial dimension, p the gas density, T the tem- 
perature, and the atomic mass was set to unity), e is the 
unit vector pointing to the position of the particle at the 
moment when it hits the sphere and 2)e is the integra- 
tion measure over angular coordinates. The measure De 
additionally depends on the relative velocity g = u — v; 
for the hard-sphere gas 2)e = ^^^ 0(g • e) de, where 6{-) 
is the Heaviside step function and de is the standard an- 
gular integration measure. The post-collision velocity v' 
of the particle can be expressed via v, e, and g: 



2e(g . 



(2) 



Equation ([I]) is applicable in the diluted limit when the 
volume fraction occupied by atoms is small: pa <C 1. 

The collision integral in Eq. (IT]) can be simplified in the 
long time limit assuming that the velocity distribution is 
isotropic and that the particle velocity greatly exceeds 
the the typical velocity of background atoms, v » ^/T. 
The last assumption will allow us to expand the collision 



integral in a power series of the small parameter y/T /v. 
It is worth noting that an initial velocity v <?C \T will 
become of order yfT after a single collision (see Eq. (l2|). 
First we treat /(v) as a function of V^ = v-v. Squaring 
(pi) we get V = V + 4(u • e)(g • e). Using this result and 
expanding /(v') = f(V') into a Taylor series up to the 
second order we obtain 



/(^') = /(^^)+4(u-e)(g.e; 
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We now insert this expression back into Eq. (IT]) and we 
are left with integrals over "Be and du. Using symmetry 
we deduce [T^] the dependence of the angular integrals 
on u and g: 



/De(u-e)(g-e)-A(u-g) 



1 



g u 
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where A, B are constants defined by integrals: 
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The integrals over u can be computed by taking into 
account that the velocity of atoms u ~ vT is asymp- 
totically negligible with respect to the particle veloc- 
ity. We thus expand the integrand into Taylor series of 
u/v. Completing the integration reduces [12] the integro- 
differential Lorentz-Boltzmann equation to the partial 
differential equation 



^=d^ 
dr dv 






T = 2a''-'ApTt 



(3) 



Interestingly, the constant B drops from the final equa- 
tion; the constant A is essentially irrelevant as it is ab- 
sorbed into the new time variable r. Repeating this pro- 
cedure it is possible to show [12] that higher terms in the 
Taylor expansion of f{V') leads to terms in Eq. (|3| that 
are suppressed by the powers of the ratio \/T/v and are 
thus negligible in the long-time limit. Equation (Js]) is 
a much simpler equation than Eq. M and can be solved 
exactly via Laplace transform [12] . In particular it is pos- 
sible to show that for a delta-function initial condition, 

f{v,T = 0) = 6{v - uo)/(^rfWo~^) '^here Qd = ff^ is 
the area of the unit sphere in d dimension, the asymptotic 
solution is [TU 
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(4) 



The same asymptotic behavior ^ holds for any initial 
condition that decays to zero exponentially or faster, oth- 
erwise the long time asymptotic behavior is still given by 
Eq. ([4| apart from the tail region which is determined by 
the initial condition [T2I. 



In two dimensions, Eqs. (|3|-(]4]) have been derived in 
Ref. [9] in the realm of a stochastic model for Fermi's 
acceleration. Even earlier, the exponential velocity dis- 
tribution was found to occur in another stochastic model 
for Fermi's acceleration [13j in which a particle is bounc- 
ing in a deforming irregular container of fixed volume; the 
velocity distribution becomes exponential independently 
of the container's shape and deformation protocol }T1] . 

We have also considered the case where the interaction 
between the particle and an atom separated by distance 
r can be described by a potential function U{r) that di- 
verges as [/ ~ e(r/ro)~'^ in the small separation limit. In 
this situation the Lorentz-Boltzmann equation reduces to 
a kinetic equation (analog to Eq. (Jsk) depending only on 
the interaction exponent A [12j : 
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where r = 2A (roei/^)"^' pTt and 7 = 2(d - 1)/A. 
Equation ^ admits a scaling solution of the form 



f 



A^(l+7)-i. (6) 



Plugging (|6| into ([5| we obtain an ordinary differential 
equation for <I>(w) which is solved to yield 



^{w) = C expl-A^if;!/^! , C 
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(7) 



Thus the asymptotic growth, {v) ~ r^, of the average 
speed and the scaled velocity distribution have universal 
behaviors that are solely determined by the interaction 
exponent A and the spatial dimensionality d. 

We now investigate the spatial behavior of the parti- 
cle which are the main results of this paper. We first 
use a heuristic argument. In the case of the hard-sphere 
interaction, the average speed is v ^ t ^ pa'^^^Tt, the 
mean- free path is ^ ~ {pa'^~^)~^ , a time interval between 
collisions is At ^ i/v, and hence the total number of 
collisions during the time interval (0, t) is 
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Using the standard random walk arguments we estimate 
the particle typical displacement 



'typ 



(.y/j^^y/Tt. 



(8) 



The striking feature is that the displacement is asymp- 
totically independent on the density of atoms and their 
size. This heuristic arguments can be extended to the 
case when the particle-atoms interaction is described by 
a potential and one finds [12_ that the displacement obeys 
the same growth law (jsl. 

To achieve a quantitative understanding of the spread 
of the particle we must analyze the joint distribution 
function /(r, v, t) since the spatial distribution function 
N{r,t) alone does not obey a closed equation. 



For concreteness, consider the one-dimensional case. 
This setting is physicaUy dubious as the particle is caged 
between two adjacent atoms, so the molecular chaos as- 
sumption (that is, the lack of correlations between pre- 
coUision velocities) underlying the Boltzmann approach 
is certainly invalid in one dimension. Nevertheless, a 
Lorentz-Boltzmann equation can be written and it ac- 
tually describes the situation when in each collision the 
scattering occurs with a certain probability (otherwise 
the particle and an atom pass through each other) [8lll5j. 
Then the joint distribution evolves according to 



df _^ df 



5/ 
dv 






(9) 



The right-hand side of Eq. ^ is the same as in Eq. ([3| 
(after using d = 1 and the proper value A — 1 in one di- 
mension), so the right-hand side is asymptotically exact 
when V ^ (u) ^ v T. Further, in this v S> \/T limit the 
convective term (wV/) can be replaced by the diffusion 
term {—DV^f). This is evident if we recall that in the 
standard Lorentz gas, the particle undergoes diffusion in 
the hydrodynamic limit (see e.g. [SHHj). In our model, 
in the large time limit the particle experiences many col- 
lisions in a short time interval when the speed remains 
almost constant. The separation between the time scale 
at which diffusion appears (few collisions) and the time 
scale at which the particle speed changes appreciably al- 
lows us to replace the convective term by the diffusion 
term. In one dimension Z? = y- [8 and thus the govern- 
ing equation becomes 



dl 
dr 



dl 
dv 






\v\ dy 
4p^T dx^ 



(10) 



where we use again r — 2pTt as the time variable. Note 
that the diffusion coefficient is proportional to velocity 
and thus on average it linearly grows in time (see Eq. Q). 
Equation (10) is asymptotically exact in the hydrody- 



namic limit, V ^ V T; the full integro-differential Boltz- 
mann equation is required for the description of the early 
stage of the time evolution when the particle has experi- 
enced only few collisions. 

In the long time limit, the joint distribution function 
f{x,v,t) should approach the scaling form 

f(x,v,t)^-^F{X,V), X^-, V=- (11) 



where x* = yTt and w* = r. The velocity and displace- 
ment reflection symmetry |12j allows us to limit ourself 
to the quadrant y > 0, X > 0. The normalization condi- 
tion for f{x, V, t) is recast in Jg°° dX J^ dV F{X, V) — I 



which explains the factor 1/4 in the scaling ansatz (11). 
It is convenient to study Eq.( |10| on the entire x-axis 
with V > 0. Performing the Laplace transform in v and 
the Fourier transform in x we find that the transformed 
joint distribution 



g{k,q,T) = 



dx( 



dve-^''f{x,v,T) (12) 




FIG. 1. Density profile for the hard sphere gas vs. the rescaled 
variable R — r/^/Tt. The numerical simulations (symbols) in 
d = 1, 2 are compared with the theoretical predictions (con- 
tinuous lines). The theoretical prediction for d = 3 is also 
shown. 



satisfies 

dg_ 
dr 



(^-Q^)| 



= -kg, 



Ap^T 



(13) 



The exact solution of this linear hyperbolic partial dif- 
ferential equation can be found using the method of char- 
acteristics to yield ^2j 



g{k,Q,T) ^ 



1 



cosh s + 4 sinh s 




Q tanh i 



4 tanh s ' 



where s = Qt and go{k,Q) — g{k,Q,T = 0). Taking 
the long-time limit (which, for rapidly decaying initial 
conditions, is equivalent to choose go=l |12| ) and explic- 
itly computing the inverse Laplace transform we obtain 
scaling function F{X, V) 



F{X,V) = 



ds 



-isX • 



-V s coth s 



sinh ; 



(14) 



We could not compute the integral ( 14 ) in a closed form, 
so we determined it numerically (Fig. pi). The spatial 
distribution {N{X) = L dVF{X, V)) admits an explicit 
expression 



N{X)^ [coshi?i 



R^ = ^X 



(15) 



The same approach holds in any dimension for a hard- 
sphere particle-atom interaction 12J. For example, in 
three dimensions 



F{R,V) 
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NiR) 



3^3 (-R3 



ds s' 
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For the hard sphere gas in one and two dimensions, 
the inhomogeneous Lorentz-Boltzmann equation (fTl) was 
simulated by stochastically updating the velocities and 
positions of 10* and 10^ particles respectively [T2]. The 
velocity distribution is in excellent agreement with the 




FIG. 2. F{X,V) for hard spheres gas in Id, Eq. (14 1. (Top 



panel) 3D-plot. (Bottom panel) Values along the lines of fixed 
V = 0.0035, 1, 2, 3, 4. The continuous lines are obtained from 
the numerical simulations while the symbols represent the 
values obtained by computing the integral (14 1. 



exponential scaling form. The density profiles are shown 
in Fig. IT] Both in one and two dimensions there is ex- 
cellent agreement with the theoretical prediction on the 
full range of the spatial coordinate. 

In conclusion we have analyzed the behavior of a very 



light particle in an equilibrium background gas. We have 
shown that in the long-time limit, the average particle 
displacement grows linearly with time and proportion- 
ally to the thermal velocity of the background atoms — 
the density of the gas, the size of atoms, and the details 
of the interaction between the particle and the atoms do 
not affect the asymptotic. The average particle veloc- 
ity also grows in a rather universal way and the scaled 
velocity distribution approaches a scaling form which is 
generically non-Gaussian (the only exception is when the 
particle-atoms interaction is described by a Maxwell po- 
tential, 7 = 1 in Eq. (Is])). For the hard-sphere particle- 
atom interaction in arbitrary dimensions, we have com- 
puted the asymptotically exact velocity distribution, po- 
sition distribution and joint velocity-position distribution 
using a combination of Fourier and Laplace transforms. 
Our analytic solutions explicitly show the lack of factor- 
ization: The joint distribution F{R, V) is not a product 
of the density N{R) and the velocity distribution ^{V). 

Our theoretical predictions are in perfect agreement 
with the numerical simulations providing strong evidence 
that our simulation scheme is correct and that the sim- 
plification of the collision integral and the replacement 
of the convective term by effective diffusion are indeed 
asymptotically exact in the limit when the particle ve- 
locity greatly exceeds the thermal velocity of atoms. 

The Lorentz model was originally suggested U as an 
idealized model of electron transport. Perhaps the most 
interesting extension of the present work is to analyze the 
quantum version of our model. 
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